Emergence of Order in Textured Patterns 
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A characterization of textured patterns, referred to as the disorder function 8{j3), is used to 
study properties of patterns generated in the Swift-Hohenberg equation (SHE). It is shown to be an 
intensive, configuration-independent measure. The evolution of random initial states under the SHE 
exhibits two stages of relaxation. The initial phase, where local striped domains emerge from a noisy 
background, is quantified by a power law decay S(/3) ~ t"' 3 . Beyond a sharp transition a slower 
power law decay of 5{j3), which corresponds to the coarsening of striped domains, is observed. The 
transition between the phases advances as the system is driven further from the onset of patterns, 
and suitable scaling of time and 5(f3) leads to the collapse of distinct curves. 

The decay of 5(f3) during the initial phase remains unchanged when nonvariational terms are 
added to the underlying equations, suggesting the possibility of observing it in experimental systems. 
In contrast, the rate of relaxation during domain coarsening increases with the coefficient of the 
nonvariational term. 



I. INTRODUCTION 

The study of spatio-temporal patterns has received 
considerable impetus from a series of elegant experiments 
and theoretical developments based on symmetry consid- 
erations. Recent experimental studies include those on 
reaction diffusion chemical systems ffl , convection in flu- 
ids Q and gases ||, ferrofluids 0, and vibrated layers 
of granular material [pi. These results have been sup- 
plemented with patterns generated in (relatively) sim- 
ple model systems The most complete theoretical 
treatments of patterns rely on the study of symmetries of 
the underlying system and those of the patterns || . Un- 
fortunately, this analysis is restricted to periodic or quasi- 
periodic patterns. A theoretical analysis of more complex 
states requires the identification of suitable "variables" 
to describe a given pattern. Examples of such variables 
include the structure factor |l(J, the correlation length 
p] |l3f and the density of topological defects G3|. In 
this paper we study properties of another characteriza- 
tion, referred to as the "disorder function" jl5],[l6| . 

The patterns studied are generated in physical systems 
(and models) whose control parameters are uniform in 
space and time; thus, they result from spontaneous sym- 
metry breaking. The simplest class of nontrivial struc- 
tures are periodic. They are typically striped, square, 
triangular or hexagonal patterns that form in perfect, 
extended arrays Ha]. To obtain periodic patterns, the ini- 
tial state of the system and/or the boundary conditions 
need to be carefully prepared. A second class consist 
of periodic patterns whose "unit cells" have additional 
structure Jl7],|l8|]. A field describing periodic arrays can 
be expanded in a few plane waves. 

The patterns described above contain a unit cell that is 
repeated on a "Bravais lattice" to generate a plane-filling 
structure. The qualitative description of the pattern in- 



volves the characterization (in terms of symmetries) of 
the unit cell and the generators of the Bravais lattice. 
For example, the unit cell of a honeycomb lattice is Dq- 
symmetric, and the Bravais lattice is generated by two 
unit vectors 120° apart. 

Quasi-periodic patterns have also been observed under 
suitable experimental conditions p9| . Their symmetries 
can be observed in Fourier space. For example, the spec- 
trum of a quasi-crystal is 10- fold symmetric j pOj . Quasi- 
periodic patterns can be described using a few "princi- 
pal" plane waves along with their nonlinear couplings. 

The bifurcations to and from a given periodic (or quasi- 
periodic) state can be studied using the "Landau equa- 
tions," which once again rely on the symmetries of the 
physical system and the pattern [ pT| . The information 
used is that, since the pattern is generated by symmetry 
breaking, a second pattern obtained under the action of 
any symmetry of the physical system has identical fea- 
tures. The imposition of this equivalence (supplemented 
by the elimination of "higher order" terms) gives the nor- 
mal form equations of the pattern. They contain infor- 
mation on aspects of dynamics of the pattern and details 
about its bifurcations 

Patterns such as those of Figs. |l]and| (which are gen- 
erated in a model system) do not belong to the classes 
discussed above. These structures, referred to as "tex- 
tured" or "natural" patterns |^2| , are observed when the 
initial states from which they evolve are not controlled. 
Similar structures are seen in small aspect ratio systems 
when the boundaries play a significant role in the creation 
of the pattern ||. 

There is no (nontrivial) global symmetry of textures; 
consequently, they cannot be characterized using sym- 
metry groups. Note also that a second realization of the 
experiment (e.g., starting from a different set of initial 
conditions) will give a pattern that is vastly different in 
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detail (such as Figs. |l|(a) and 0(b)). In spite of these dif- 
ferences, one can clearly recognize similarities between 
distinct patterns. For example, the correlation length 
and the density of topological defects of the two tex- 
tures shown in Fig. |l| are similar. In contrast patterns 
generated under other external conditions (e.g., Fig. ||) 
have different characteristics. A theoretical treatment of 
textured patterns requires a "configuration independent" 
description. 

FIG. 1. Two patterns generated by evolving random initial 
states via the Swift-Hohenberg equation for 1600 time units. 
The parameters used were D = 0.1, e = 0.2, v = 2 and fco = 1. 
The initial states consisted of white noise whose intensity var- 
ied between and 10~ 3 . Periodic boundary conditions were 
imposed on the square domain of 256 x 256 lattice points, the 
length of whose sides are (487r/fco). 

In Section II, we introduce the disorder function 8(0) 
whose definition was motivated in part by the argu- 
ment leading to the derivation of Landau equations; for 
patterns generated in uniform, extended systems, 8(0) 
is required to be invariant under rigid motions of a 
labyrinthine pattern jn^pjj . In section III, we briefly 
describe the method for evaluating the disorder function 
from (typically noisy) grid-values of the field. It relies 
on a method to approximate a continuous function from 
values given on a grid referred to as the method of "Dis- 
tributed Approximating Functionals" (DAFs) 0| . 

The main results of the paper, which include properties 
of the disorder function, and its application to provide a 
quantitative description of the relaxation of the patterns 
from an initially random state are presented in Sections 
IV and V. The underlying spatio-temporal dynamics is 
given by the Swift-Hohenberg equation (SHE) and one 
of its variants. We first provide evidence to support the 
claim that the disorder function consists of intensive, con- 
figuration independent variables. The use of 8(0) shows 
that pattern relaxation occurs in two distinct phases sep- 
arated by a sharp transition. We also study changes in 
the relaxation profile when the system is driven further 
away from the onset of patterns. In Section V, we discuss 
the effects of adding nonvariational terms to the SHE. 

FIG. 2. Two patterns generated by evolving a random ini- 
tial state via the Swift-Hohenberg equation for 2400 time 
units. The parameters used for the integration were D = 0.01, 
e = 0.4, v = 2 and ko = 1/3. The initial states consisted 
of white noise whose intensity varied between ±10 -2 . The 
length of each side of the square is (487r/feo). 



II. THE DISORDER FUNCTION 

Textured patterns observed in experimental systems 
and those shown in Figs, [l] and |^ can be described 
by a scalar field v(x) which is smooth, except perhaps at 



the defect cores. However, unless the patterns are triv- 
ial (e.g., perfect stripes, target patterns) the analytical 
form of the field is unknown. Consequently, it is diffi- 
cult to determine a set of "configuration independent" 
characteristics of structures generated under similar con- 
ditions. We instead impose a weaker requirement, that 
the characterizations remain invariant under the action 
of the symmetries of the underlying physical system; i.e., 
translations, rotations and reflections p5| . Rather sur- 
prisingly, the measures so defined have similar values for 
distinct patterns such as those shown in Fig. [l| 

The most significant feature of labyrinthine patterns is 
that they are locally striped; in a suitable neighborhood 
w(x) ~ sm(k-x), where the modulus ko(= |k|) of the 
wave vector does not vary significantly over the pattern. 
Structures generated in experiments and model systems 
include higher harmonics due to the presence of nonlin- 
earities in the underlying system; they only contribute to 
the shape of the cross section of stripes. In order to use 
the simplest characterization of textures, we eliminate 
the second and higher order harmonics by the use of a 
suitable window function in Fourier space. For experi- 
mental patterns (which do not have periodic boundary 
conditions) this is a nontrivial task, and a method to 
implement it is described in Ref. |24j |. 

The simplest local field that is derived from w(x) and 
whose value remains the same under all rigid motions is 
its Laplacian Ati(x). Terms such as A"w m (x), though 
invariant, are difficult to extract from an incompletely 
sampled field (typically given on a square lattice). The 
requirement that perfect stripes be assigned a null mea- 
sure (they are not disordered), coupled with the local 
sinusoidal form of the (filtered) pattern implies that the 
lowest-order field relevant for our purpose is (A+fcg)i>(x). 
The family of measures, referred to as the disorder func- 
tion, is defined by 



8(0) = (2-0 



Jda\(A + kZ)v(x)f 
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where < |w(x)| > denotes the mean of |t>(x)|, and 8(0) 
has been normalized so that the "intensive variables" 
8(0) = 8(0)/ J da are scale invariant. The moment is 
restricted to lie between and 2 for reasons discussed be- 
low. Local deviations of the patterns from stripes (due 
to curvature of the contour lines [^5)) contribute to 8(0 
through the Laplacian, while variations of the width of 
the stripes contribute via the choice of a "global" ko. 

8(0) depends on the choice of the wave-numer ko of 
the "basic" stripes. Analysis of striped patterns u s t (x) = 
Asin(\i ■ x) and target patterns Wf(x) = Acos(kr) show 
that 8(1) is minimized when ko = k = \k\. Studies of tex- 
tured patterns from model equations indicate (see Fig. |^) 
the presence of a unique minimum of 8(1). We use the 
minimization of 8(1) as the criterion for the choice of the 
wave- number fcg in Eqn. ([!]). For patterns generated us- 
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ing the Swift- Hohenberg equation g], fco is very close to 
the wave-number obtained by minimizing the "energy" 
p5| . We also find that our estimation of fco is far more 
robust (i.e., smaller variation between distinct patterns) 
than that evaluated from the power spectrum. This is 
presumably because wavelength variations and curvature 
of contour lines at each location of the pattern contribute 
to the computation of 5(0). 

The variation in k over the pattern (traditionally de- 
fined to be the half width of the Structure Factor 27 1) 
can be estimated using the variation of 6(1) with the 
wave- number fc (see Fig. [|) . In the remainder of the pa- 
per we define Ak to be the distance between fc-values for 
which 5(1) is twice the minimum value p6[ . Analysis of 
textures shows that Afc is a configuration- independent, 
intensive variable. 



domain wall contains curvature of the contour lines and 
variations of the stripe width; consequently it has non- 
zero disorder. 5(0 for a single domain wall is a mono- 
tonically increasing function of the angle 9 between the 
stripes of the two domains plfl . Thus 5(0) provides in- 
formation absent in characterizations such as the corre- 
lation length. The disorder function for a target pattern 
u(x) = acos(kor) is known |l5| ], and is used to determine 
the accuracy of the numerical algorithms. For target pat- 
terns, the integral in the numerator diverges as (2 — /3) _1 , 
and leads to limiting the range of (3(< 2), and to the in- 
troduction of the prefactor in the definition of 5(0. 

The weights of distinct characteristics of a texture 
(e.g., domain walls, defects, variations of the stripe- 
width, etc.) depend on the moment (3. In particular, 
the contribution to 5(0 from a domain wall vanishes as 
(3 — ► 2, and the limit is proportional to the number of 
targets in a pattern. The effects of noise on the calcula- 
tions are minimal. For example, addition of 10% white 
noise typically changes 5(0 by less than 2%. 
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FIG. 3. The behavior of 5(1) as a function of the 
wave- number k for a labyrinthine pattern, fco is estimated 
to be the (unique) minimum of the curve. The width Afc is 
defined to be the distance between the fc-values at which 5(1) 
reaches twice its minimum. 

Observe that our choice of fco is arbitrary in one sense; 
we could have chosen to minimize 5(0) for any fixed j3 
to determine fco. However, the observed variations in fco 
are insignificant. Alternatively, we could have evaluated 
each 5(0) by minimizing it with respect to fc. We choose 
not to implement this scheme because of the need to es- 
timate only one free parameter (i.e., the wave- number) 
for a given pattern. 

For a perfect set of stripes the function 5(0 = 0. A 



III. DISTRIBUTED APPROXIMATING 
FUNCTIONALS 

The critical requirement for a good estimate of 5(0 
is a sufficiently accurate determination of the Laplacian. 
Calculation of derivatives of a field from values given on 
a lattice is a delicate task, especially in the presence of 
noise. A technique, utilizing what are referred to as Dis- 
tributed Approximating Functional (DAFs), has been 
introduced recently, to fit analytically or approximate a 
continuous function from known values on a discrete grid 
p3[ . Unlike typical finite difference schemes, it estimates 
the function and its derivatives using a range of neighbor- 
ing points (~ 40 in our case); consequently, the required 
computations are much less sensitive to noise. 

The most useful for our application has been a class 
of DAFs for which the order of accuracy of the fit is 
the same both on and off the grid. (This is in contrast 
to interpolation, which forces the fit to be exact on the 
grid, but always leads to intertwining about the func- 
tion off the grid, thus leading to less accurate estimation 
of derivatives.) Alternatively, the method is designed so 
that there are no special points. If the labels on the grid 
points are erased after DAF-fitting of a function, it be- 
comes essentially impossible to identify the points that 
were on the grid. The most general derivation of the 
DAFs is via a variational principle [p9| , yielding 

9daf(x) = ^ I( x : x k )g(x k ), (2) 

k 

where k labels the grid points, g(x k ) are the known input 
values of the function (which may contain noise), and the 
sum is over x k such that \x k ~ x\j 'A < i?, A being the 
lattice spacing. For suitable I(x,Xk), the function and 
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its derivatives are evaluated to a comparable accuracy 
p9| . This proves crucial in the evaluation of the disorder 
function. 

The calculations presented are carried out using the 
"Hermite DAF" @, defined by 

Ae-* 2 ^ 2 / iy' i 
I(x,x k ) = I(x - x k ) = = V - 7 -H 2j (z), 

V^TT j=0 \ 4/ J! 

(3) 

where z = (x — £fc)/<7v2 and H n (z) is the n tfl Hermite 
polynomial. The Gaussian weight (of width a) in Eqn. 
(||) makes I(xi — Xk) highly banded, reducing the com- 
putational cost of applying the DAF to data. The DAF 
representation of derivatives of a function known only on 
a grid is given by 

(0)^W=E^^)^), (4) 

which can be evaluated either on or off the grid. In the 
continuum limit, the derivative of the DAF equals (ex- 
actly) the DAF of the derivative |(J. 

The DAF approximation to a function that is sampled 
on a square grid (x m , y n ) can be obtained using the two- 
dimensional extension 

I((x,y), (x m ,Vn)) = Ix(x,x m )I Y (y,y n ) (5) 

of the approximating kernel |5l| . Thus to estimate (say) 
, Eqn. (0) needs to be applied in the y-direction (with 
/ = 0) and along the 2-direction (with 1 = 2). (The 
application of the DAF operators in the two directions 
commute and can be carried out in any order.) 

As M — > oo, I(x.Xk) — > 5(x — Xk), and the DAF ap- 
proximation qdaf{x) becomes exact. With finite M, the 
sum on the right side of Eqn. (3) becomes a polyno- 
mial (of order M), resulting in gDAF (x) being smooth. 
Note that gDAF(x) can be considered to be a weighted 
running average of the signal. The Gaussian width a 
determines the effect of neighboring points in the DAF- 
approximation. The range R is chosen to be sufficiently 
large that the terms ignored in Eqn. (2) are negligi- 
ble. (The Gaussian weight included in the definition of 
I(x, Xk) guarantees the decay of these terms.) 

IV. PATTERNS GENERATED USING THE 
SWIFT-HOHENBERG EQUATION 

The patterns analyzed in the paper are obtained from 
periodic fields w(x, t) generated by integrating random 
initial states through a modified Swift-Hohenberg equa- 
tion (SHE) H|| 

d t u = D(e-(ko + A) 2 ) u - 7 w 3 - v(\7uf. (6) 



The parameters D, kg, and 7 can be eliminated through 
suitable rescaling of t, x, and u respectively, e measures 
the distance from the onset of patterns. The results for 
the variational case [y = 0) are presented in this Section 
and those for the nonvariational case (y ^ 0) are given 
in the next. 

The initial fields for the integration were random num- 
bers in a predetermined range. The time evolution is im- 
plemented using the Alternating Direction Implicit algo- 
rithm jj4|. Each nonlinear term iV[w(x, t)] is expanded 
to first order in 5u = u(x, t+St)—u(x., t), thus linearizing 
the equations in u(x.,t + 8t). Updating the field involves 
the inversion of a penta-diagonal matrix. The typical 
time step used for the integration, At ~ 0.1, was chosen 
so that the higher order terms in Su are insignificant. We 
have confirmed the robustness of the integration by com- 
paring (in a few cases) the results with those done for a 
smaller time-step (At ~ 0.001). 



A. Properties of the disorder function 
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FIG. 4. The curves 5(0) for patterns generated at two dif- 
ferent sets of control parameters. The lower bunch consists of 
curves for four patterns at the first set of control parameters 
(Fig. 1) while the upper bunch consist of those for a second 
set of control parameters (Fig. 2). (b) shows the same plots 
with a logarithmic vertical scale. 
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The form of the measures 5(f3) were deduced by re- 
quiring invariance under rigid motions of a single pat- 
tern. Are these limited restrictions sufficient to yield 
characterizations that can delineate the observed "com- 
monality" in distinct patterns generated under identical 
conditions? Surprisingly, it appears to be the case. Fig. 
U shows the disorder functions for several patterns. The 
curves bunched at the bottom show 8(0) for four tex- 
tures (two of which are shown in Fig. generated at 
fixed control parameters. 5 ((3) appears to have captured 
the commonality of these distinct patterns. Structures 
generated in the Gray-Scott model |3^] and in a vibrated 
layer of granular material 0] exhibit similar properties 




The next question is if 8(0) can differentiate between 
patterns with different visual characteristics. Fig. ^ 
shows two structures obtained from the SHE for a sec- 
ond set of control parameters. They have characteristics 
that differ from patterns of Fig. Q e.g., they contain 
smaller domains and a larger density of defects. 5(0) 
for four such textures are bunched together on the upper 
curves in Fig. ^. The significant separation of the two 
sets of curves (e.g., the values of 5(1) between the two 
sets is about 25 times larger than the average difference 
between curves within a set) confirms the ability of 8(0) 
to quantify the differences of the two groups of patterns. 



FIG. 5. The values of 5(0.5), 5(1.0) and 5(1.9) for periodic 
patterns generated in domains of different sizes. The areas 
of the domains are 367r x 367r, 367T x 72n, 72n x 72tv, and 
144-7T x 144-7T. In each case random initial states are evolved 
for 8000 time units under the SHE with control parameters 
given in Fig. 2. The results indicate that 8(0) are intensive 
variables for labyrinthine patterns. 

The disorder function quantifies the characteristics of a 
labyrinthine pattern using the local curvature of the con- 
tour lines and the wavelength variations, which typically 
increase with the (visual) disorder of a texture. Thus, 
8(0) is able to quantify the observation that patterns of 
Fig. H are more disordered that those of Fig. [|. 

Next, we provide evidence to substantiate the claim 
that 8(0) are intensive variables for labyrinthine patterns 
such as those shown in Figs. Q and ^ |3jJ. This is done 
by comparing values of 8(0) for patterns (with periodic 
boundary conditions) of several sizes. The sizes of the 
domains chosen are 367r x 367r, 36tt x 72ir, 72tt x 72tt, 
and 144-7T x 1447T, and each pattern is generated by inte- 
grating a random initial state (with amplitude between 
±1CT 2 ) for a time T = 8000 under the SHE. The results, 
shown in Figure 0, give the mean of 10 patterns for each 
domain size (except the largest where only 5 patterns 
were used). The results indicate that 5(0.5), 5(1.0), and 
5(1.9) are intensive variables, and the corresponding 8(0) 
are extensive variables. 



B. Relaxation of patterns 
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The characterization of textures using 8(0) finds one 
useful application in the study of the relaxation from 
an initially random state. Fig. |^ shows several snap- 
shots of a relaxing pattern. During an initial period 
(t < T ~ 800) the local domains emerge out of the 
random background and the mean intensity < |w(x, t)\ > 
nearly reaches its final value. The subsequent evolution 
due to domain coarsening is very slow. These qualitative 
features are repeated in multiple runs under the same 
control parameters. 

Figure shows the behavior of Afc, 5(0.5), 5(1.0) and 
5(1.9) for the evolution shown in Fig. EL The curves re- 
main identical (except for small statistical fluctuations) 
for different realizations of the experiment; i.e., the disor- 
der function captures configuration independent aspects 
of the organization of patterns. The relaxation clearly 
consists of two stages, with a sharp transition in 5(0) at 
t = T Q 0JTJ. 
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FIG. 6. Several snap-shots of the relaxation of random initial state whose intensity is between ±10 -2 under the SHE with 
D = 0.01, e = 0.4, v = 2 and fco = 1/3. An initial phase (t < 800) when the local striped patterns are being formed is followed 
by domain coarsening. 
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FIG. 7. The behavior of (a) Afc, (b) 5(0.5) (c) 5(1) and (d) 5(1.9) during the evolution shown in Fig. 6. During the 
initial phase 8((3) decays (approximately) like ' . During the second phase the scaling is nontrivial; e.g., 5(0.5) ~ t~ 09 , 
5(1.0) ~ t -015 and 5(1.9) ~ t -019 . The transition between the two phases occurs around t = 800. 



During the initial phase, the time evolution of 5(1) 
changes smoothly from a logarithmic decay to a power 
law 5(1) ~ t~ 71 , where 71 w 0.5. Corresponding de- 



cay has been observed in the width of the structure factor 
p0| . The scaling is "trivial" in the sense that for other 
"moments" S((3) - @. The decay of 5(0) appears 
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to be associated with the L ~ 1 5 growth of domains in 
non-conserved systems [ [37] . 

The second phase of the relaxation (due to domain 
coarsening) exhibits a more complex behavior. The mo- 
ments <5(0.5), 6(1) and <5(1.9) behave approximately as 
i~ - 09 , t~ - 15 and i~°- 20 respectively, indicating the pres- 
ence of "non-trivial" scaling Since the relative con- 
tribution of isolated defects increases with (3 (Section II), 
the slower decay of 5(1.9) (compared to 6(1) 19 ) suggests 
that changes in the density of defects is less significant 
than the reduction of curvature of the contour lines [[L0| . 
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ing e. Similar results are observed for all values of the 
moments 0. 
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FIG. 9. The scaling function for the relaxation obtained by 
the rescaling t' = te and 5' = 5e~^ . The data correspond to 
e values of 0.01, 0.05, 0.1 0.4 and 0.8. 

Suitable scaling of variables, including t — > et can be 
used to eliminate e from the SHE. Hence we expect that 
the rescaling t — ► et and 6(1) — > e~?6(l) will lead to col- 
lapse of the curves shown in Fig. ||. This is indeed the 
case as seen from the scaling function (Fig. ph. 



FIG. 8. The evolution of S(l) for patterns generated by 
integrating random initial states under the SHE for several 
values of e. Each curve is an average of five runs. For clar- 
ity, the standard deviations are shown only for one parameter 
value. For distinct e, 5(1) exhibits identical behavior during 
initial growth of domains, and decays at the same rate dur- 
ing the coarsening phase. The other moments 5(f3) exhibit 
similar behavior. 



C. Changes in the relaxation with e 

Figure ^| shows the behavior of 6(1) during the relax- 
ation of random initial states under the SHE for several 
values of e, all other parameters being fixed. The initial 
decay of 6(1) and the rate of decay during the second 
phase are seen to be independent of e. Furthermore, the 
transition between the two phases advances with increas- 



V. RELAXATION IN NONVARIATIONAL 
SYSTEMS 

In this Section we discuss properties of 5(f3) when the 
spatio-temporal dynamics is nonvariational. The absence 
of an underlying "energy" of the dynamics suggests a 
faster relaxation, since the system cannot be constrained 
by "metastable states" during the evolution. The behav- 
ior of the disorder function confirms this expectation. 

Figure [l(] shows the behavior of 6(1) for the organiza- 
tion of a random field under a nonvariational SHE (i.e., 
!/^0). The decay of 6(1) remains the same (as the anal- 
ogous variational dynamics) during the initial relaxation 
(Fig. 0) and becomes faster during the coarsening phase. 
As the coefficient v of the "nonvariational term" in Eqn. 

increases (the value of the remaining coefficients re- 
maining the same) , so does the relaxation rate during the 
coarsening phase (Fig. 0). 
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The wave- number (obtained by minimizing 5(1) over 
ko) relaxes to a value (ko = 0.61) that is larger than the 
corresponding one for the variational case (ko = 0.59). 
Such a deviation was observed earlier in Ref. [T^ ] where 
it was suggested that ko (for the nonvariational case) cor- 
responds to the zero-climb velocity of isolated dislocation 
defects. 
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FIG. 10. The behavior of 5(1) for the evolution of a ran- 
dom initial state under the nonvariational modification of the 
SHE. The parameters of the SHE are the same as given in 
Fig. 6 except for v — 2.0. The initial decay is identical to the 
variational case while the coarsening phase exhibits a faster 
decay of 5(1). 



local striped domains emerge out of the noisy background 
and their amplitudes saturate close to their final value. 
This behavior is described by a logarithmic decay fol- 
lowed by a power law decay of the disorder 5(1) ~ t~? . 
The scaling is "trivial" in the sense that the decay of the 
remaining moments satisfy 5(0) ~ 5(1)^. The second 
phase of the relaxation corresponds to domain coarsen- 
ing and is a much slower process. The scaling during this 
phase is nontrivial. 

As the system is driven further from the onset of pat- 
terns (as measured by the parameter e) the duration of 
the initial phase is reduced. However, the rates of decay 
of the disorder function for the two phases remain un- 
changed and rescaling of time by e and of 5(0) by e~2^ 
leads to a scaling collapse. 
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VI. DISCUSSION 

We have used the disorder function 5(0) to character- 
ize properties of textured patterns and their relaxation 
from initially random states. The disorder function was 
defined by requiring its invariance under rigid motions of 
a single texture. It was found to be identical for multiple 
patterns generated under similar external conditions; i.e., 
5(0) is configuration independent. We provided evidence 
to confirm that the moments are intensive variables. In 
addition, the disorder function can differentiate between 
patterns with distinct characteristics. 

The evolution of initially random states under the 
Swift-Hohenberg equation is conveniently described us- 
ing 5(0). The relaxation consists of two distinct stages 
separated by a sharp transition. During the initial phase, 



FIG. 11. The behavior of 5(1) for several values v. Ob- 
serve that the decay during the growth phase remains the 
same, while the decay during the coarsening phase increases 
with v. 

The addition of nonvariational terms to the spatio- 
temporal dynamics leads to several interesting observa- 
tions. The decay of disorder during the initial phase is 
unchanged, and appears to be a model independent fea- 
ture. Thus, one may expect to observe it during relax- 
ation of patterns in experimental systems. The expecta- 
tion of a faster relaxation in nonvariational systems (due 
to the absence of "potential minima") is seen only dur- 
ing domain coarsening. This rate of relaxation is system 
dependent and increases as the coefficient of the nonva- 
riational term. 
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There is very little theoretical understanding of the ob- 
served behavior of the disorder function. The decay of 
S(f3) during the initial phase of pattern relaxation is rem- 
iniscent of analogous behavior in the XY- model p8| , |39"|] , 
and appears to correspond to the L ~ 1 2 growth of do- 
mains in nonconserved systems [ p7| . 

We would like to thank K. Bassler, D. K. Hoffman, R. 
E. Jones, D. J. Kouri and H. L. Swinney for many stim- 
ulating discussions. This work was partially funded by 
the Office of Naval Research and the Energy Laboratory 
at the University of Houston. 
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